Deposit Growth in the Wetting of an Angular Region with Uniform Evaporation 
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Solvent loss due to evaporation in a drying drop can drive capillary flows and solute migration. 
The flow is controlled by the evaporation profile and the geometry of the drop. We predict the 
flow and solute migration near a sharp corner of the perimeter under the conditions of uniform 
evaporation. This extends the study of Ref. 6, which considered a singular evaporation profile, 
characteristic of a dry surrounding surface. We find the rate of the deposit growth along contact 
lines in early and intermediate time regimes. Compared to the dry-surface evaporation profile of 
Ref. 6, uniform evaporation yields more singular deposition in the early time regime, and nearly 
uniform deposition profile is obtained for a wide range of opening angle in the intermediate time 
regime. Uniform evaporation also shows a more pronounced contrast between acute opening angles 
and obtuse opening angles. 

PACS numbers: 47.55.Dz, 68.03.Fg, 81.15.-z 



I. INTRODUCTION 



Evaporative contact line deposition, the "coffee-drop 
effect" has been the subject of several recent papers 
[1EIESHISH0- The physical problem originates 
from a simple phenomenon of everyday life: when a drop 
containing a solute such as coffee dries on a surface, the 
solute is driven to the contact line, forming a characteris- 
tic ring pattern. This simple phenomenon is potentially 
important in many areas of both scientific and industrial 
applications [3, llO, UjL LL2II ■ This evaporation mechanism 
can create very fine lines of deposition in a robust way 
that requires no explicit forming. Further, it is a way of 
concentrating material strongly in a quantitatively pre- 
dictable way. Lastly, it creates capillary flow patterns 
that can be useful for processing of polyatomic solutes 
like DNA 0El[lil. 

One striking aspect of this deposition phenomenon is 
its dependence on the shape of the droplet. This de- 
pendence was recently explored by Popov and Witten 
[fjLl, who studied corner-shaped drops. Here the liquid 
region on the surface has the form of a sector of arbi- 
trary opening angle a (Fig. Such shapes contrast 
strongly with the circular drops treated in previous stud- 
ies II II II 11 Popov and Witten found that this 
difference in shape led to striking differences in evapo- 
rative flow and deposition near the apex of the drop. 
Both the flow and the deposition profile showed singular 
power-law behavior as a function of distance r from the 
tip. These power laws vary in a predicted way with the 
opening angle, but arc otherwise universal. The growth 
of the deposition also shows several different predicted 
behaviors in three defined time regimes. Thus by manip- 
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ulating the shape of a droplet one has extensive control 
over its deposition properties. 

In this paper we further explore the range of control 
possible in evaporative deposition. This study is moti- 
vated by a finding from Ref. 6: the power laws govern- 
ing the deposition depend on the evaporation conditions 
in the vicinity of the drop. Ref. 6 considered the usual 
diffusion-controlled evaporation conditions, in which the 
liquid is surrounded by a dry surface and the evapora- 
tion rate is limited by the diffusion of vapor away from 
the drop. In these conditions concentration n(r) in the 
air obeys Laplace's equation with n at the surface set by 
the fixed saturation concentration. As in the analogous 
electrostatic problem, the gradient of n diverges at the 
edge and at the tip. The resulting divergent evaporation 
profiles contribute strongly to the controllable deposition 
properties found in Ref. 6. 

One may readily alter these evaporation conditions, 
and strong differences in the deposition are expected to 
result. We consider a condition that contrasts strongly 
with the singular evaporation treated in Ref. 6, viz. uni- 
form evaporation. This contrast is illustrated in Fig. ^ 
Uniform evaporation can be created by surrounding the 
drop by a wet surface instead of a dry one, as illustrated 
in Figure IB. The wet-surface evaporation case allows 
us to discern the role of the evaporation profile in pro- 
ducing the power-law behaviors revealed by Ref. 6. The 
wet-surface evaporation case is also mathematically sim- 
pler to treat than the dry-surface case of Ref. 6. This 
allows us to make cleaner predictions with fewer approx- 
imations for the wet-surface case. 

Our work is to a large extent a mere application of 
the general theory worked out in Ref. 6. We shall follow 
the same approach and exposition for the current case 
that was used in Ref. 6. Before working out the quanti- 
tative behavior we describe the experiment qualitatively. 
A droplet of linear dimension L is forced to have a cor- 
ner shape over part of its perimeter. One may fix the 
edge of the drop e.g. by making shallow scratches in 
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FIG. 1: Sketch of the two different evaporation conditions 
contrasted in the text. Both pictures show a small liquid 
drop of size L and maximum thickness h m with part of its 
edge constrained to an angular corner shape by means of e.g. 
two scratches on the surface. Figure A shows a drop on a dry 
surface. Magnified tip region indicates the profile of evaporat- 
ing flux across the surface by a row of arrows. Shading indi- 
cates the wet region. The flux diverges at the edges. Figure B 
shows a similar drop surrounded by a wet surface. Magnified 
region shows the uniform evaporating flux. 



the surface. These serve to pin the contact line at the 
scratch. We choose L smaller than a few millimeters so 
that gravitational effects on the droplet 's shape are mi- 
nor. We denote the maximal thickness of the drop by 
h m and keep the droplet volume small enough to assure 
that h m -C L. Once the evaporation starts, the volume 
diminishes at a constant rate, and thus h m decreases lin- 
early with time. At some final time tf this thickness 
extrapolates to zero. Like Reference 6, we restrict our 
attention to time much less than this t f , and it will show 
later that for t Ci/, in the early drying stages, the time 
dependence of L and h m can be actually ignored. Our 
predictive power is strongest for this regime. 

The shape of the droplet during this thinning process is 
dictated by its surface tension. Near the edges, the local 



reduction in volume owing to thinning is much smaller 
than the local loss of volume to evaporation. Thus a 
flow towards the edge is needed in order to replace the 
evaporative loss. Any solute suspended in the fluid is 
carried along by this flow. The asymptotic flow near the 
tip is minimally influenced by the bulk of the drop, and 
thus this flow can be readily calculated. Knowing this 
flow profile, one may deduce how much solute should be 
carried to a given point on the edge in a given time t. 
This amount grows in a characteristic way with time and 
with distance from the tip. Our aim is to see how the 
deposition is influenced by opening angle and how much 
this deposition differs from the dry-surface results of Ref. 
6. The main difference is in the opposite direction from 
what one might expect. We find that the wet-surface 
evaporation leads to deposition that is more concentrated 
at the tip than the dry-surface deposition of Ref. 6. This 
is despite the diverging evaporation at the edges and tip 
produced in the dry-surface case. 

The paper is organized in parallel with Ref. 6. First, 
we review the basic physical model and its mathematical 
framework, which was introduced in Ref. 6. Next, the 
system is solved both analytically and numerically: the 
flow field is described, and its asymptotic properties stud- 
ied. We obtain the power law of the deposition rate in 
early time regime and intermediate time regime. Then we 
compare our results with those of the dry-surface case. 
More discussion and conclusions follow in the last sec- 
tions. 



II. PRINCIPLES GOVERNING THE 
DEPOSITION 

Following Ref. 5, we consider a droplet of very dilute 
suspension bounded in an angular region of opening angle 
a, as illustrated in Fig.^ We use cylindrical coordinates 
(r, ip, z) defined in Fig. [21 with the range < r < oo 
and —a/2 < ip < a/2, to describe our system, z is the 
coordinate normal to the substrate. We first present the 
governing equations in generality and then specialize to 
the regime of interest: thin drops undergoing slow flows. 

The equilibrium shape of the surface is dictated by the 
minimization of surface energy. The minimum energy 
surface has uniform mean curvature H. Specifically, if 
the surface tension of the liquid is a and the pressure 
difference across the surface is Ap, then the balance of 
normal forces dictates: 

Ap = -2Ha. (1) 

In the sequel, we shall represent the surface by its thick- 
ness profile h(r, ip, t). The mean curvature H depends 
on the local derivatives of h, to be specified below. 

As indicated in the Introduction, evaporation induces 
a flow towards the edge of the drop. Denoting the local 
velocity by u(r, ip, z, t) with in-plane component u s , it 
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FIG. 2: Sketch of the angular drop showing geometrical quan- 
tities used in the text 



is useful to define a depth-averaged velocity field v 



=- / Usdz. 



(2) 



The condition for local mass conservation may be stated 
in terms of this v: 



p 



(Vh) 2 + d t h = 0. 



(3) 



Here p is the density of the fluid, and Jo the mass loss 
per unit projected area and time at the point in question. 
In Ref. 6, this Jo is a strong function of position, and it 
diverges at the contact line. In the present work, it is a 
mere constant. By itself, this condition is not sufficient 
to determine v. Ultimately, v is determined by New- 
ton's equations on each fluid element. We shall consider 
the creeping flow regime where forces are in near equilib- 
rium and acceleration plays a negligible role in Newton's 
equations. Then Newton's equations reduce to the Stokes 
equation: 



Vp = ryV u, 



(4) 



where p is the fluid pressure, and rj is the dynamic viscos- 
ity. The lack of inertia implies that v is a potential-like 
flow, as shown below. 

Further physical considerations and simplifications are 
needed to solve the system analytically. Firstly, as we 
are considering a thin drop, there is a separation of the 
vertical and horizontal scales in this problem, and several 
simplifications follow 6] . The pressure inside the drop p 



does not depend on the z coordinate d z p = 0. The surface 
of the drop should have a small slope \Vh\ <C 1. And 
the z-derivatives of flow dominate, i.e., \d z Ui\ ^ \d s Ui\, 
where s represents any direction parallel to the substrate 
plane, and Ui refers to any velocity component. Under 
these considerations, (@J has the form |(J 



0, d z u s 



V s p = r]d zz u s 

With boundary conditions: u s | s:= o 
obtain 

Vp (z 2 
u s = — I — — hz 

v V 2 

With expressions J2J an d ©i we have |(j 
v = -—Vp. 



(5) 
0, we 



(6) 



(7) 



Thus v/h 2 can be represented as a gradient of the scalar 
potential as announced above. The flow in the drop 
clearly depends on the relative importance of surface 
forces and viscous forces. This importance is ordinar- 
ily characterized by the capillary number Ca — vrj/a. 
For water-like fluids, the capillary number is small when- 
ever v <C 100 m / s. Ordinary evaporating flows are 
in the range of 10 -5 m / s; accordingly, viscous forces 
may be considered as weak in comparison to surface ten- 
sion. Thus we anticipate that the shape of the drop 
is nearly the equilibrium shape in the absence of the 
flow. To establish this formally and systematically, we 
express the pressure and height as expansions in the cap- 
illary number: p = po + (Ca)pi + (Ca) 2 pi + ... and 
h = ho + (Ca)hi + (Ca) 2 /i2 + By using these expan- 
sions in Eqs. QJ, © and J7J, we find the lowest order 
results El 
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V-(AgVV) = -— -d t h 



v =h Q Vtp, 



(8) 



(9) 



(10) 



where ip — — pi(Ca)/ (3rj), and the leading term po does 
not vary with local coordinates and is only a function 
of time t. Thus, one can use Eq. (JHJ to determine the 
equilibrium drop surface shape ho, then solve Eq. @ with 
respect to tf)(r,ip,t), and finally determine the velocity 
field from Eq. (|1U|) . To simplify the notations, we will 
write po — p a tm as Ap, and ho as h in the rest of the 
paper. 

Eq. (JSJ introduces a length scale into the problem, i.e., 
the mean radius of curvature R(t) = a/Ap. Since the 
evaporation rate is constant in time, the droplet volume 
and thus its thickness decrease linearly in time until the 
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terminal time tf. This means @ that the mean curvature 
R(t) oc L 2 /h m can be written as 



R(t) = 



FL 



1-t/t/ 



(11) 



where Ri is the initial mean radius of curvature, and tf 
is the total drying time. Thus in the early drying stages 
(t -C tf), which is the only case we are to consider in this 
paper, the time dependence of R can be ignored, and the 
shape of the drop can be assumed to vary with time adi- 
abatically (i.e., slowly compared to all other processes). 
We may generally treat R(t) as constant Ri in the rest 
of the paper, keeping in mind its time dependence (|1 1|) . 



III. BEHAVIOR OF THE DROP NEAR THE TIP 
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If written explicitly in terms of h(r, tp, t), equation (JSJ 
reduces to the Poisson equation under the assumption of 
small slope of the surface: 



(12) 



with boundary conditions h(r,—a/2) = h(r,a/2) = 
h(0,tp) = 5]. Here a/Ap = R(t) w R t is the mean 
radius of curvature, and we are only interested in the 
asymptotic limit r <C R(t). This problem is quite classi- 
cal |l7j . however, its results are somewhat unexpected [5J. 
The leading term of the solution in the limit r <C R(t) is 
different for acute and obtuse angles and can be obtained 
by either solving the small-slope (horizontal) equation 
(|12|l or by the series expansion of the full equation 
By both methods, the lowest order term in r of the series 
expansion of the surface shape h was found in Ref. 5: 



h(r,cp,t) 



R(ty 



where 



- , , If cos 2ip 

h{<p) = - £ - 1 

4 V cos a 

hUp) — C (a) cos — , 
a 



0<a< -, 



7T 7T 

-, 77 < " < 7r - 
a 2 



(13) 

(14) 
(15) 



Exponent v as a function of the opening angle a is shown 
in Fig. 03 The acute-angle solution {T3J, is self- 

similar and independent of the remote boundary con- 
dition specifying the surface shape in the bulk of the 
drop. The obtuse-angle solution (|13fl . (|15(l has a pre- 
factor C (a) , which does depend on the remote boundary 
condition in the bulk of the drop 01 an d has the form 
0: 



C(a) 



1 



4a - 2ir 



C + O(a 



(16) 



FIG. 3: Relation between the exponent v and the opening 
angle a 



where Co is independent of a and is determined by the 
remote boundary condition. Clearly, the leading-order 
solution 113|) , 114|) , and (|15fl diverges when a approaches 
the right angle from either side. 

The notable difference between acute and obtuse open- 
ing angles in the main order results (|T3)l . (CH|) . and i(T5|) 
was hrst discussed in Ref. 17. For acute opening angles, 
a local similarity solution is valid, i.e., the solution that 
has no arbitrary coefficients and is independent of the 
remote boundary conditions in the bulk of the drop that 
need to be further prescribed. Physically, the surface 
shape of the drop near the tip is controlled locally and 
is independent of the physical conditions far from the 
tip. For obtuse opening angles, however, the local simi- 
larity does not hold, and the surface shape is no longer 
controlled locally. Accordingly, the coefficient C (a) de- 
pends on the remote boundary conditions. This contrast 
between acute and obtuse opening angles leads to differ- 
ent flow behavior and different deposition properties as 
will be shown later. 

The dependence of C (a) on the remote boundary con- 
dition can be illustrated in the following way. C (a) 
can be expressed in terms of the linear dimension L 
and the maximal thickness h m , which characterize the 
global shape of the drop. Along the bisector ip = 0, 
let r — > L and h —> h m , and using approximation 
R(t) ~ Ri oc L 2 /h m , one can obtain from Eqs. i|13fl and 
131: 



do ) V ( A 



(17) 



Clearly, C (a) is small since v < 2 for obtuse opening an- 
gles, and its time dependence is weak in the early drying 
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stages and can be ignored. 

For numerical purposes, Co in (|16fl will be set to unity, 
as was done in Ref. 6. This is justified by the fact that 
when a approaches tt/2 from above, the divergent term 
l/(4a — 2tt) dominates, and the properties of the sys- 
tem become independent of Cq. In the opposite limit, 
a = 7T, the tip of the angular region can be chosen at 
any point on the contact line due to the symmetry, and 
physical properties of the drop are again independent of 
Cq. Although for arbitrary obtuse opening angle Cq is 
controlled by the remote boundary conditions, numerical 
studies have shown consistent results that did not vary 
substantially with Cq. 

The special case a — tt/2 invites further explanation. 
According to expressions (|14|) and l|15fl . function h(tp) di- 
verges when a approaches tt/2 from either side. This di- 
vergence is artificial, however [{| . This issue is explained 
in some detail in Appendix B. In the rest of the work, 
we will treat a = tt/2 as the limiting case and will keep 
in mind the possible divergence of expressions l|14|) and 

03. 
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B. Reduced pressure 

Now we are ready to solve Eq. We are interested 
only in the asymptotic behavior when r — » (or r <C R), 
and in this limit dth oc r v as we know from Eq. (|13fl . 
Thus, dth can be safely dropped compared to the con- 
stant term Jo/ p on the right-hand side of Eq. @ • Phys- 
ically, in this asymptotic limit during the early drying 
stages the fluid mass transport due to the gradient of the 
flow flux is uniquely balanced by the evaporation from 
the surface locally at each moment, with the mass change 
brought about by the local height change being a higher 
order small quantity that can be ignored. 

The asymptotic solution for the reduced pressure 
ip(r, ip, t) of Eq. © can be expressed as: 



ip(r, (p, t) = 



Jo 



(18) 



p m 3 - 

where the exponent of r is determined by simply counting 
the powers of r on the left side of Eq. @. From Eqs. © 
and (|18|) . we find explicitly the differential equation that 
■0 (ip) should satisfy Q 



d 2 ip 3dhd^ ~ 
T~ 2 + ~ ~T~ ~~j 2(3u — 2)ip 



1 

h 3 



T7T- (19) 



dip 2 ' h dip dip 

This equation depends implicitly on the opening angle, 
as —a/2 < ip < a/2, with v and h depending on a as 
shown in Eqs. (|T4"|) and ljl"5j) . 

The boundary conditions associated with Eq. I|19|) need 
to be clarified. Firstly, we expect the flow field to be 
totally symmetrical with respect to the bisector ip = 0, 
and therefore tp should be even in ip: 



dip_ 

dp 



= 0. 



(20) 



FIG. 4: A, B: The reduced pressure function ^(ip) for repre- 
sentative acute and obtuse opening angles. A: a = 7r/4, B: 
a — 37r/4; C, D: The regularized function x(p), defined in 
Appendix A as x(v) ~ h 2, ij){ip), corresponding to A and B 



Secondly, the outer boundary condition at ip — a/2 can 
be identified explicitly as 



'3 dtp 
dp 



= 0. 



(21) 



tfi=a/2 



In Appendix A, we justify this boundary condition, and 
discuss how Eq. (fTT)|l may be regularized. 

The boundary problem of Eq. (|19|) with boundary con- 
ditions 1)20(1 and (|21() is complete and has a unique solu- 
tion. We show the numerical solutions ip (ip) for typical 
opening angles a — it/ A and 37r/4 in Fig.^J The analyti- 
cal solution to this boundary problem and its asymptotics 
are also discussed in Appendix A. 

Eq. (|19fl is solvable analytically for special opening an- 
gle a = 7r. If written in terms of the Cartesian coordi- 
nates (where the x axis is the contact line and the y axis 
is the bisector), the reduced pressure function ip can be 
found from Eqs. (|15fl . (|18|) . and (|19[) . and has the form 



iP 



Jo 1 1 

P C 3 (tt) y 



(22) 



tp=0 



Since is does not depend on x, this result is fully consis- 
tent with the symmetry of the system when a = n. 

For < a < ir, no special angle exists to reduce the 
complexity of the equation. Angle n is the only opening 
angle where we can obtain analytic results in a closed 
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form. For a = tt the simplicity is well anticipated, since 
this limiting case has no apex at all. Without an apex, 
all points on the boundary are equivalent, and most of 
the resulting properties follow by symmetry. We will use 
the exact solution at the opening angle tt to test our 
numerical results and analytical asymptotics. 



IV. RESULTS FOR THE PHYSICAL 
PROPERTIES OF THE SYSTEM 



A. Flow field 



In polar coordinates, the velocity field (|10[) can be ex- 
pressed explicitly: 



v =v r f + v v ip, 



-(3,-2)^ T 



Jo ( r 
p \R(t) 



p \R(t) 

dp 



(23) 

(24) 
(25) 



As we assume that solute particles carried by the fluid 
move with the same velocity as the flow itself (this as- 
sumption was actually confirmed both theoretically and 
experimentally, and a theoretical estimate can be found 
in [lg), the trajectory of each particle is identified as 
the streamline of the flow field, and can be obtained by 
integrating the velocity field [6j: 



dr v r 
rdip Vtf 



dtp 
dtp 



(26) 



If the particle eventually arrives at (ro, a/2) on the con- 
tact line, the trajectory reads as |(j 



f a/2 - ( d$ 
r(p) = r exp ( (3v - 2) J 



dp' . (27) 



In Fig. [3 we show the flow field configuration for the 
opening angle a = tt/2 computed numerically with the 
acute- angle expression (THE)! . 

These flow trajectories are independent of the dimen- 
sional quantities like the evaporation rate J, and they 
depend only on the opening angle @- Thus, the trajec- 
tory of the moving solute particle does not depend on 
how fast the drop evaporates and how thin the liquid 
layer is. 

We now consider the asymptotic properties of the 
streamlines. For p — > a/2, velocity component v r 
vanishes due to the outer boundary condition, which 
amounts to the statement that no mass flows inward or 
outward along the contact line, while component v v re- 
mains finite due to the asymptotic behavior of dip /dp 




FIG. 5: Streamline configuration for the opening angle n/2 



(All), and goes as 
Jo 

P 



R(t) 



dh /a\ 

'dip' 12 J 



(28) 



or, if written explicitly by employing results (|14|) and 

(113, 



2Jo 1 / 



P 
Jo 



tana \R(t) 
a / r 



p irC(a) \R{t) 



1 — TT I a 



0<a<|, 



— < a < 7r. 
2 



(29) 



In particular, v v seemingly vanishes as a — > 7r/2. How- 
ever, in fact it does not. We explore this point in Ap- 
pendix B. Consequently, as p — > a/2, dr/dp — > 0, and 
the streamlines are perpendicular to the contact line. 

In the opposite limit, p — > 0, the streamlines reach far 
into the bulk of the drop, and function r (p) is divergent. 1 
Behavior of the streamlines depends on the asymptotic 
properties of ip near the bisector, which, according to Eq. 
(|19fl and boundary condition l|2U|). can be written as 



i> - </J(0) 



2(3^-2)^(0) 



1 



(30) 



The boundary value tp(0) requires the complete solution 
of Eq. I|19fl . The divergent part of the integration in 127|) 
is then uniquely determined by the asymptotic form l|30() 
since (Eq. (50) in Ref. 6): 




dp' 



In 



a 



2(3v-2)-n 2 2p 



where 



h 3 (0)ip(0) 



(31) 



(32) 



1 Since v/h 2 is a potential flow, the streamlines r{tp) may end only 
where the reduced pressure if) assumes local extreme values, i.e., 
at the contact line and at infinity (in the bulk of the drop). 
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Therefore the streamlines in this limit scale as 

r^r Q {£)\ ^0, (33) 

with 

(34) 



2(3^-2)-k 2 ' 

and ro = r(a/2) is the contact line distance. 

For special opening angle tt, it is straightforward to 
obtain the streamline equation and the velocity compo- 
nents in terms of Cartesian coordinates from Eqs. (|22|l . 
(|21|), and (J2SJ: 



x = r 



0. 



1 Jo 
'C(tt) p ■ 



(35) 



(36) 



Accordingly, we have k 2 (tt) = 1 and e(7r) = 1. The 
streamline configuration and velocity are fully consistent 
with the symmetry of the system, since the position of 
the apex is no longer well defined when a = tt. 

We do not have exact analytic results for exponents 
k 2 (a) and e (a) for arbitrary opening angle a, since equa- 
tion (|19|l can not be solved in closed form to obtain ip(0). 
One has to solve the boundary value problem i|A13|) nu- 
merically to fix ^(O), and determine k 2 and e in terms 
of the opening angle a. We show the dependence of k 2 
and e on the opening angle a in Figs. |5] and [7| for com- 
parison, we also include results found in Ref. 6 for the 
case of dry-surface evaporation. Numerical results are in 
agreement with the analytic result we have found for the 
special case a = tt. 

The exponent e determines the asymptotic behavior of 
the streamlines near the bisector ip = 0. According to 
Eq. 133|) . the distance between a streamline and the bi- 
sector scales with ip as tpr(ip) oc t^ 1 " 6 in this limit. As 
is apparent in Fig. \7\ in the case of wet-surface evapo- 
ration e is equal to unity for a — 7r/2 and tt, and hence 
tpr (ip) remains constant asymptotically. Geometrically 
this means that streamlines run parallel to the bisector 
when ip — > 0. This result also follows directly from our 
analytic solution (|35|l for a = tt and from Fig. [5] for 
a = 7r/2. For a = tt this geometric property of stream- 
lines is well anticipated from the symmetry of the sys- 
tem; for a = tt/2, however, it is not obvious. For acute 
opening angles, we find e < 1, and therefore the asymp- 
totic distance decreases when ip — > 0, and the streamlines 
converge toward the bisector as r — > oo. The incoming 
particles are moving along the trajectory away from the 
bisector. For obtuse opening angles, we have e > 1, and 
therefore the asymptotic distance increases when <p — > 0, 
and the streamlines diverge away from the bisector as 
r — > oo. Now the incoming particles first move along the 
trajectory toward the bisector, reach a minimal distance, 
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FIG. 6: Dependence of parameter a 2 on opening angle a. The 
solid line corresponds to the case of uniform evaporation. The 
dotted line, obtained in Ref. 6, corresponds to the dry-surface 
evaporation with quadratic profile (Eq. (25) in Ref. 6). 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



20 40 60 80 100 120 140 160 180 



a (degrees) 



FIG. 7: Dependence of the stream line asymptotic exponent 
e (when ip — > 0) on the opening angle a. The solid line cor- 
responds to the uniform evaporation case. The dotted line, 
obtained in Ref. 6, corresponds to the dry-surface evaporation 
with quadratic profile (Eq. (25) in Ref. 6). 



8 



and then turn away toward the contact line. For dry- 
surface evaporation [(|, streamlines always diverge away 
from the bisector, since e > 1 in this case for all opening 
angles (Fig.[7J). 



B. Solute transfer and deposit growth 

As in the dry-surface evaporation case (Eqs. (55) and 
(56) in Ref. 6), one can now calculate the time it takes 
for the solute particles initially located at (r, tp) (where 
< (fi < a/ 2) to move along the streamline to the contact 
line (r a/2). We use Eqs. ((23) and J23 to write: 



t = 



a/2 

tp v f> 



rd(p 



(37) 



W2 exp L(3^-2)/; /2 ^)(fr 
'"/ — ^ 



where 



to 



Jo R{ty- 



j rt 



(38) 



is a combination of the system parameters that has a 
dimension of time. In the early stages of the drying pro- 
cess, we can approximate R(t) by and to is therefore 
independent of time. For each streamline indexed by 
contact line distance ro, and for each time t, there exists 
a unique ip{ro,t) determined by Eq. H37fl . such that all 
the solute located in the area bounded by neighboring 
streamlines indexed by r and r + dr in the domain 
fiTo^t) < C < a /2 reaches the contact line and becomes 
part of the deposit within time t. Let us denote the mass 
accumulated between r and r + dr at the contact line 
by time t as dm(ro,t). 

To study solute transfer and deposit growth, we want 
to understand the deposit distribution along the contact 
line, as well as its growth rate with time. We can consider 
the amount of solute d(dm) /dt arriving at the contact line 
during time dt through the flow tube bounded by neigh- 
boring streamlines indexed by ro and ro + dro, as shown 
in Figure|S] (compared to Ref. 6, where a global approach 
was employed, and the integral quantity with respect to 
time dm was computed). If we assume that the initial 
concentration of the solute is constant Co everywhere in 
the drop, then d(dm)/dt can be expressed as: 



d(dm) 
dt 



(r ,t) = c h(r(r ,ip(r ,t))) \v(r) x dr\ . (39) 



Again, we can start with the special opening angle 7r, 
for which analytical results can be obtained straightfor- 
wardly. In Cartesian coordinates, the surface shape i|15|) 
reads as h = C(7r)y, the velocity is given by Eq. (|36|l . 
and y = v y t. Expression (|39|l can be simplified as 



d{dm) 
dt 



(r ,t) = c h(f(r Ql ip(r ,t))) \v y \dr , (40) 




To+C/Ab 



FIG. 8: Qualitative sketch of the local approach. The so- 
lute in the shaded area will be arriving at the contact line 
through the flowtube at time t, and this corresponds to 
d (dm/dt) (rot). 



from where a power law can be found: 



d 
dt 



dm 
dr 



Co 



C(tt) 



(41) 



The deposition rate for a = n does not depend on ro, 
which is anticipated, since the position of the apex is no 
longer well defined. Thus, for a = tt there is a unique 
power law in the whole domain of tp and for all times t 
in the early drying stages. 

For arbitrary opening angles a other than 7r, however, 
the deposition rate can not be calculated analytically 
without knowledge of the closed form of ipi'fi)- Instead, 
we have to analyze ()37|l as well as other relevant quanti- 
ties asymptotically in two different limiting cases: <p — > 
and if — > a/2. These two different asymptotic regions 
correspond to two different time regimes, when the de- 
position rate follows different power laws, as was first 
introduced in Ref. 6: 

early time regime: t <C to, when only solute par- 
ticles initially located near the boundary can reach the 
contact line and become part of the deposit. Proper- 
ties of deposition in the time regime are governed by the 
asymptotic ip — > a/2. 

intermediate time regime: to <C t tf, when so- 
lute particles initially located near the bisector are able 
to reach the contact line. This time regime is governed 
by the limit tp — > 0. The condition f < t/, where tf is 
the total drying time, means that we are still considering 
early enough drying stages, where our model applies. 

As argued in Ref. 6, the separation of two time regimes 
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works worse when the opening angle a increases toward 
7r. This can be readily seen in our case: as shown in (|41ll . 
the deposition rate follows the same power law through- 
out the early drying stages when a = 7r, and the two time 
regimes are indistinguishable. 



1. Deposit growth in the early time regime 

In the early time regime, only the solute particles ini- 
tially located near the contact line contribute to the mass 
deposition. 

As shown earlier, the velocity component v r vanishes 
in the limit ip — > a/2, and streamlines are perpendicular 
to the contact line. Expression (|39[1 has a very simple 
form: 



d(dm) 
dt 



oc v v hdro, 



a 

^ 2- 



(42) 



In this limit, the asymptotic forms are h oc rg(a/2 — ip) 
(Eq. (H3), v v oc r\~ v (Eq. (gSJ), and t(ip) oc r%(a/2-(p) 
(Eq. and therefore Eq. l|4^1) can be written in the 

form of a power law: 



d 
dt 



dm 
dr 



oc tr. 
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where 



(3 = 1- v. 



(43) 



(44) 



One can also use the approach employed in Ref. 6 to 
find the power law of the deposition rate (compared to 
Eq. (60) in Ref. 6): 



dm 
dr 



2J3 



(r .t) oc t r„ 



(45) 



where (3 is again given by (|44|l . which is in agreement 
with result (|43() and the result we found for special angle 
7r (|4ip. The dependence of the exponent (3 on the opening 
angle is shown in Fig. where the exponent obtained 
in Ref. 6 for the case of dry-surface evaporation is also 
included for comparison. 

The relation l|44|l can be understood in the following 
way. Let S be the distance from the contact line. When 
d is small, mass conservation demands hv oc JS, and 
hence v oc (dh/dS)^ 1 , i.e., near the contact line the 
velocity should be inversely proportional to the slope 
of the drop surface \dh/d6\. According to Q13[l. near 
the contact line h = r^h(ip) « r v h(a/2 — S/r), and 



\dh/d5\ 



(dh/d(p) 



is proportional to r 1 
and Eq. (glj follows. 



. Therefore, the velocity 
so is the mass deposition rate, 



2. Deposit growth in the intermediate time regime 

In the intermediate time regime, we need again to find 
the power law of d(dm)/dt. By analyzing expressions 
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FIG. 9: Dependence of the exponent /3 of the contact line 
distance ro on the opening angle a in the power law dm/dro 
(Eq. 1451 1 in the early time regime. The solid line corresponds 
to the uniform evaporation, and the dotted line corresponds 
to the dry-surface evaporation with the quadratic profile (Eq. 
(25) in Ref. 6). 



(E3, HI, (123, 03, and in the limit tp -> 0, we 
find that relevant physical quantities assume the follow- 
ing asymptotic forms: 



r oc r$Lp 



(46) 



v r oc -r^tp-^ 1 -^, v v oc rl- v y l -< l - v \ (47) 



{«r V'. 



(48) 



Then, according to Eq. ll3l?|l , 

d{dm 



dt 



oc h \v(r) x df\ — h \v v dr — v r rd(p\ . (49) 



By assuming dt = along the direction of df, we can de- 
rive the relation dip oc — (tp/r^) dr^ from Eq. I|48|l . Using 
Eqs. 1461) and (|47|l . expression (|49|l can be simplified as: 



d ( dm 



dt \ dro 1 

Together with H48JI . we finally obtain: 

d ( dm\ s _ 1 

cct rj, 



dt \dro 



(50) 



(51) 
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FIG. 10: Dependence of the exponent 5 of time t on the open- 
ing angle a in the power law dm/dro (Eq. 1541 ) in the in- 
termediate time regime. The solid line corresponds to the 
uniform evaporation, and the dotted line corresponds to the 
dry-surface evaporation with the quadratic profile (Eq. (25) 
in Ref. 6). 



FIG. 11: Dependence of the exponent 7 of the contact line 
distance ro on the opening angle a in the power law dm/dro 
(Eq. JSU) in the intermediate time regime. The solid line 
corresponds to the uniform evaporation, and the dotted line 
corresponds to the dry-surface evaporation with the quadratic 
profile (Eq. (25) in Ref. 6). 



where 



5 = 1 



u(3u-2) 



(52) 



and 



7=1 



3^-2 



(53) 



In this limit, we can also follow Ref. 6 to compute the 
power law of dm/dro, and we find the deposition rate to 
be (compared to Eq. (63) in Ref. 6): 



dm , . I 
dm 



(54) 



V. COMPARISON WITH THE DRY-SURFACE 
EVAPORATION CASE 

Popov and Witten [|| considered a general evaporation 
rate J(r, ip) of the asymptotic form 



J 



-x 



(55) 



The wet-surface case considered here corresponds to /x — 
1 and A = 0. The results of the previous section may be 
obtained by setting fi = 1 and A = in the general ex- 
pressions of Ref. 6. In some cases this limit permits sim- 
pler expressions and more explicit solutions as we have 
seen. For the dry-surface evaporation case studied in 
Ref. 6, A = 1/2 and fi is an explicit function of the open- 
ing angle a (Fig. 1 12(1 reflecting the singular behavior of 
the Laplacian vapor concentration field |lSj . Compared 
to the dry-surface evaporation, the uniform evaporation 
rate yields different deposition properties in both early 
and intermediate time regimes. 



with the same scaling exponents given by l)52|) and (|53ll . 
Again Eqs. I|51l) and l|54l) are in agreement with each 
other and exponents l|52|l and (|53|) are in agreement with 
the values of exponents we obtained for special angle ir. 
We show the dependence of the exponents 5 and 7 on the 
opening angle in Figs. ITUl and llll where we also include 
results of Ref. 6 for comparison. 



A. Early time regime 

In this regime, only the asymptotic form (|55|l of the 
evaporation rate J when ip approaches ±a/2 matters. 
While /x is relevant to the dependence of the reduced pres- 
sure function ip on coordinate r, the exponent A governs 



11 




0.0-1 , 1 1 , 1 , , , 1 

20 40 60 80 100 120 140 160 180 

a (degrees) 

FIG. 12: Relation between the exponent fj, and the opening 
angle a, in the case of diffusion-controlled evaporation from 
an angular region 



more fluid and mass to the sides and thus away from the 
apex. 2 

Furthermore, the 1/tq dependence of the deposition 
rate in the early time regime for acute opening angles 
seems to violate the integrability and suggests a logarith- 
mic divergence: an arbitrarily large fraction of the mass 
may accumulate within an arbitrarily small distance from 
the apex. To resolve this possible singularity, we note 
that l/r dependence occurs not at all times, but only 
at early times, which are almost never observed near the 
vertex. Indeed, for each time t one can define a crossover 
length 

r*(t)= (tRr'jY". (56) 

According to the definition of the time scale to <|38fl . the 
early time condition t <C to corresponds to the condition 
ro 3> r* (t) in terms of the location of that regime along 
the contact line. Thus, at each time t, the early time 
regime is observed only away from the vertex, and the 
areas near the vertex are always in the intermediate time 
regime. The natural cutoff r*(t) actually saves the mass 
deposition at the tip from being logarithmically infinite. 



the singular property of the reduced pressure function 
near the contact line (a direct observation is that, due 
to nonzero value of A = 1/2 in the dry-surface evapora- 
tion, the order of divergence of ip is higher by 1/2 in that 
case). For the deposition rate in the early time regime, 
the power law found in the dry-surface evaporation case 
0, as well as our result l|45[) . is uniquely determined by 
the local singular properties of the reduced pressure func- 
tion at the contact line, and those power law exponents 
have simple algebraic expressions in terms of /x, A, and 
v. 

Both results of Ref. 6 and our results show that the 
exponent of time t in the power law of deposition rate 
is independent of the opening angle a in the early time 
regime, although in our case the deposition process goes 
slower with exponent 2 instead of 4/3 in Ref. 6. As re- 
gards the dependence on the contact line distance ro, 
both results of Ref. 6 and our results (Fig. |UJl show that 
f3 always remains between —1 and for obtuse open- 
ing angles a, and the integrability of the singularity at 
r = holds. As argued in Ref. 6, despite having larger 
deposition rate, the vertex of the sector does not dom- 
inate the sides, and the deposition accumulation at the 
vertex is not qualitatively different from the deposition 
accumulation on the sides for obtuse opening angles. 

For acute opening angles a, although (3 remains be- 
tween — 1 and for the dry-surface evaporation, it is con- 
stantly — 1 in our case. It seems that we have a more con- 
centrated deposition pattern with a nonsingular form of 
the evaporation rate. The singular flow of the dry-surface 
evaporation rate J may deflect the streamlines towards 
the contact line more than in our case, and hence drive 



B. Intermediate time regime 

Distinctive power laws characterizing this regime (such 
as the configuration of the streamlines and the deposition 
rate) depend on the entire range of <p, including the limit 
ip — > 0, which dominates properties of these power laws. 
An important parameter in this regime is the value of the 
reduced pressure function at the bisector ^(0), or the pa- 
rameter k 2 (Eq. (|32|l ). Near the bisector, both the evap- 
oration rate J and function ip are regular, and in order to 
find ip(0) one needs to solve the differential equation l|19|) 
in the full domain of ip from to a/2. Therefore, in the 
intermediate time regime relevant quantities depend on 
the functional form of the evaporation rate J and surface 
shape h in the entire domain of ip and require solution 
of the main equation (|19fl . in contrast to the early time 
regime where only the singular behavior at the contact 
line matters. 

The most interesting result in this regime is certainly 
the qualitatively different behavior of physical quantities 
in the geometry of the acute opening angles versus ob- 
tuse opening angles. Our numerical result for k 2 (Fig. 
EJl is not qualitatively different from that found in Ref. 
6, as k 2 changes dramatically at a = n/2 in both cases. 



2 We actually plotted the streamlines for both dry-surface and 
wet-surface evaporation cases, however, we did not notice any 
substantial qualitative difference. This result is important and 
should be provable experimentally, since one should be able to 
measure mass accumulation at the sides and at the vertex and 
compare the two evaporation cases. 
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This is understandable, because the non-smoothness at 
a = 7r/2 is uniquely due to the crossover of the expo- 
nent v in the leading order term determining the surface 
shape (Eqs. i|14|) and (JT5J ) , while the evaporation rate J 
behaves smoothly in both cases. However, changes for 
other quantities that take place at the right opening an- 
gle are qualitative in the wet-surface case, unlike in the 
dry-surface case. 

The flow field configuration near the bisector as <p — > 0, 
is governed by the exponent e shown in Fig. [7| In the 
case of dry-surface evaporation, e is larger than 1, except 
for a — 7r, and the streamlines asymptotically diverge 
from the bisector. In contrast, our wet-surface results 
show that the streamlines asymptotically converge to the 
bisector (e < 1) for acute opening angles, while for obtuse 
opening angles they still diverge. In particular, for the 
critical angle ir/2 streamlines run parallel to the bisector 
line, which is not anticipated intuitively. 

For the power law of the deposition rate in the inter- 
mediate time regime, we obtained the asymptotic form 
(|51Jl . The result for the exponent S (a) (Fig. fol- 
lows the same qualitative pattern as the result in Ref. 6, 
although it is everywhere larger, so that the deposition 
process goes slower, like in the early time regime. 

For dry-surface evaporation, 7(a) remains negative, 
monotonically increases with a, and becomes at a = tt 
(compatible with the symmetry). In the wet-surface 
evaporation case, 7(a) follows a richer pattern. Phys- 
ically, as 7 is strictly larger than —1, the integrability 
property holds. However, while 7 remains negative for 
obtuse opening angles, it becomes positive for acute an- 
gles, and therefore the deposition accumulation is in favor 
of the sides rather than the vertex as in Ref. 6. In par- 
ticular, as 7 goes to at a = tt/2 and its absolute value 
remains small nearby, it seems that the wet-surface evap- 
oration case yields a relatively uniform deposition pattern 
for a wide range of opening angles in the intermediate 
time regime. 

The exponent 7 controlling mass deposition is closely 
related to the trajectory exponent e. From expressions 
(|r?4|l and 1)53(1 one can obtain a simple relation: 



(57) 



1 



7 



which is also implicit in Ref. 6. Intuitively, in the in- 
termediate time regime, virtually all the solute between 
the bisector and the contact line is swept to the contact 
line and becomes part of the deposit, and therefore 7, 
the exponent of the contact line distance ro indexing the 
streamlines, should be related to the geometric distribu- 
tion of the streamlines near the bisector away from the 
vertex. For a = n and 7r/2, exponent e = 1, and the 
streamlines are uniformly distributed near the bisector. 
Therefore, the solute is uniformly carried to the contact 
line by the flow, and the deposition rate should be in- 
dependent of the contact line distance, hence 7 = 0. 
For other angles, the streamlines become unevenly dis- 
tributed away from the vertex. When e > 1 , they diverge 



away from the bisector as ip — * 0, the smaller contact line 
distance tq is, the nearer the corresponding streamline to 
the bisector, and more solute will be carried to the spot 
along the streamline. Therefore the deposition is in favor 
of the vertex, and 7 is negative. When e < 1, 7 is positive 
by the same argument. 

To be more precise mathematically, in the intermedi- 
ate time regime to ^ t the deposition is uniquely deter- 
mined by the streamline configuration near the bisector 
line (this is in contrast to the early time regime, where the 
deposition rate is closely related to the slope of the sur- 
face shape near the contact line). The amount of solute 
deposited near ro is controlled by the width of the gap be- 
tween adjacent streamlines indexed by ro and r + Ar 
near the bisector line, which is proportional to rAip in 
the limit ip — » 0. Near the bisector line, the streamline 
indexed by ro can be expressed as r oc rQ(p~ e (Eq. i|46|) ). 
If we consider a small patch r = const near the bisector, 
and study the intersections of those streamlines with this 
patch, we have 



= dr cx tp £ G?ro — erop c 1 dip, 



and therefore 



dip <p 
dr er 



(58) 



(59) 



According to the above argument, Eqs. (|46(l and l|59|) 
yield 



dm dtp 1 1 / 1 
—. — oc — — oc -ro 
dro dro e 



(60) 



and the relation i|57|) follows immediately by comparison 
of the last expression with Eq. I(54(l . 



VI. DISCUSSION 

The theoretical framework, first established in Ref. 6 
and studied here in continuation, captures the essential 
mechanism of the deposit growth, but does not take into 
account a number of additional effects that can modify 
the deposition, and some restrictions and shortcomings 
remain [|J. 

The crossover at the opening angle a = n/2 (as mani- 
fested in the surface shape h (I13|) . I|14l) and (|15fl ). which 
is certainly independent of the evaporation rate, is quite 
a subtle point of the theory. As shown in Appendix B, 
in principle, modifications to the results obtained in this 
work as well as in Ref. 6 by the asymptotic analysis are 
needed in the neighborhood of the right opening angle 
[n/2-Aa,n/2 + Aa], with Aa ~ (tt/4) \\n(r / R t )\^ (Eq. 
(|B9|0 . The thinner and flatter the drop is, and the closer 
to the apex, the better our results apply. 

Further, we believe the power law exponents obtained 
by the asymptotic analysis are exact, except for a possi- 
ble logarithmic modification at a — tt/2, where crossover 
of exponents could happen. Exact properties of the other 
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experimentally testable physical quantities (the veloc- 
ity components, the reduced pressure function tjj) in the 
neighborhood [tt/2 — Act, tt/2 + Aa] could be in princi- 
ple interpolated, as shown in Appendix B, since all the 
physical properties of the system should depend on the 
opening angle a continuously. 

Another interesting observation is the form of the evap- 
oration rate. In this paper, we considered the uniform 
evaporation, which is simpler than the dry-surface evap- 
oration studied in Ref. 6. The form of the evaporation 
rate is not arbitrary, and apart from other physical re- 
strictions, it should be compatible with the symmetry of 
the system. One such consideration would be the follow- 
ing: in the limit a — > tt the position of the apex is no 
longer well defined, and the deposit rate should not de- 
pend on the contact line distance r$ in both early and 
intermediate time regimes; therefore (3(tt) = 7(71") = 0. 
Combined with the general expression for ft (Eq. (61) in 
Ref. 6), this condition demands 

A(7t)+m(7t) = 1, (61) 

which is satisfied in both cases: 3 for the dry-surface evap- 
oration A = 1/2 and /i(7r) = 1/2, and for the uniform 
evaporation A = and \i = 1. Although 7(71) = 
in our case, this did not hold exactly for the exponent 
7 in the dry-surface case, since the approximation of J 
used for the numerical solution in Ref. 6 was not in exact 
conformity with the symmetry. In contrast to f3, where 
(3(tt) = can be traced back to the simple relation (|6H . 
the exponent 7 is related to the parameter k 2 , which de- 
pends on the solution of the main equation i|19|) in the 
entire domain of tp. It is interesting to see how the differ- 
ential equation (|19f) . which is merely a statement of local 
conservation of mass, together with the proper form of 
the evaporation rate, yields the deposition properties in 
full conformity with the symmetry of the system. 

In general, various properties of the evaporating drop 
in the intermediate time regime depend on the mathe- 
matical structure over the whole domain of ip, and more 
thorough analytical treatment of the main equation 119|) 
is certainly appealing. The criticality of the opening an- 
gle a — tt/2 demands extra attention: Why do some 
physical properties differ between acute angle and ob- 
tuse angle? Why was not this separation so apparent in 
Ref. 6, where the same h entered? How is this separa- 
tion related to the form of the evaporation rate? Naively, 
these questions can be readily addressed by saying that 
the criticality of the right angle is uniquely due to the 
crossover of the leading terms in the full expansion of 



The physical origin of this condition is very simple: for the open- 
ing angle a = 7T, the evaporation rate should depend only on the 
normal distance to the contact line rsin(a/2 — which re- 

duces to r(a/2 — \ip\) when \ip\ — ► a/2. Thus, the exponents of 
r and (a/2 — \<p\) must be equal for this opening angle, and the 
general expression 1551 immediately yields \x — 1 = —A. 



the surface shape h at a = tt/2. When the evaporation 
rate is uniform, it does not introduce any further singu- 
larity, and thus helps to retain the trace of the criticality 
of h in the resulting physical quantities and phenomena. 
In the case of the dry-surface evaporation, the stronger 
dependence of J on coordinates may overshadow the co- 
ordinate dependence of h, and the latter may appear less 
significant in the results. In particular, it seems that the 
change of the exponent \i in the evaporation rate J is 
in such a direction as to compensate for the change of 
exponent v in the surface height h. More mathemati- 
cally rigorous treatment needs to be done to make this 
argument clearer. 

Experimentally, it is interesting to note some possible 
applications of our results. Our work shows that wet- 
surface evaporation at early times with acute opening an- 
gle achieves the greatest concentration towards the apex. 
Accordingly, one could actually achieve a great concen- 
tration of mass by allowing the evaporation to occur for 
a short time, then allowing the dissolved solute to diffuse 
and equilibrate, then allowing another bit of evaporation, 
and so forth. In this way, one could approach the behav- 
ior of having a finite fraction of the mass within some 
small distance of the apex. Another interesting aspect 
is the nearly perfect uniformity of the deposition for the 
intermediate time regime with opening angle a ~ 7r/2. 
This kind of uniform deposition may be useful, especially 
when a small amount of the concentrated substance is 
sufficient. For example, a dilute solution of reagents can 
be concentrated strongly at the contact line, thereby in- 
ducing a chemical reaction there. The evaporation mech- 
anism assures that the concentration is a known function 
of the position and the initial dilution. Likewise, trace 
amounts of solute can be rendered more easily detectable 
by causing them to concentrate at a contact line. 



VII. CONCLUSION 

The wet-surface evaporation of an angular drop yields 
surprisingly rich and potentially useful behavior. This 
behavior complements the previously studied work on 
dry-surface evaporation [(|. Though our case lacks the 
distinctive singular evaporation of the dry-surface evapo- 
ration case, remarkably, it leads to a stronger focusing of 
solute towards the apex. Further, it can create two qual- 
itatively different types of flow, according to whether the 
opening angle is acute or obtuse. The deposition profile is 
remarkably uniform for the intermediate times when the 
opening angle is close to a right angle. Now that these de- 
position properties have been established, they may well 
prove useful. For example, they provide a means of con- 
centrating trace solutes in a liquid in a rapid and quan- 
titatively predictable way. They also create distinctive 
capillary flow fields and distinctive concentration profiles 
of solute. We expect this kind of microscopic, singular, 
evaporative flow to play an increasing role in the tech- 
nology of small scale material synthesis, processing and 
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analysis. 
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APPENDIX A: REDUCED PRESSURE 
EQUATION: OUTER BOUNDARY CONDITION 
AND REGULARIZATION 

In this appendix, we provide a justification of the 
outer boundary condition i|21|) . and we also show how 
the boundary problem can be regularized. In Ref. 6, Eq. 
I(19jl was solved using the global mass conservation con- 
dition 



r a/2 

pi \v r \ hrdip 
'-a/2 



a/2 



Jnrdrdip, 



(Al) 



J-a/2 



where v r = h 2 d r ip is the radial component of the flow 
velocity 6J. However, Eq. i|19|) in itself is an expression of 
local mass conservation following Eq. J3J). It is therefore 
interesting to find out explicitly how Eq. IjAljl represents 
new information that can restrict the solution. Here we 
explain the origin of the new information and show that 
it takes the form of an explicit boundary condition. 

With the expressions of h and ip, Eq. (|A1|) can be sim- 
plified as 



a/2 



(A2) 



The local mass conservation implicitly stated in Eq. 119(1 
allows us to express ip in terms of the derivative dip / dip: 



2(3i 



2)fh 3 = h 3 pt 
dtp* 



- 2 dh dip 



dp> dp 



(A3) 



Then integration by parts, together with the boundary 
condition converts the integral condition (| A2|) into 
the boundary condition l|21l) . 

We now show that condition 121(1 does in fact restrict 
the general form of ip and fixes the order of its divergence 
near the contact line. To demonstrate this, we study the 
asymptotic behavior of ip. As both h and tp are even 
in (p, we can consider only the domain < (p < a/2. 
We restrict our attention to the limit tp — > a/2 and keep 
in mind that h oc (a/2 — ip) in this limit. From the 
expressions of h (|14f) and (|15fl . as well as Eq. (|19|) . we see 
that tp necessarily diverges in this limit, and the third 
term on the left side of Eq. (|19|) . which is of the lowest 
order of divergence, can be neglected in the asymptotic 
analysis. Thus, Eq. 1(19(1 reduces to: 



d 2 ^ 3 dh dtp 1 
dip 2 h dip dip h? 



0. 



(A4) 



Eq. (I!) can be solved analytically. First, consider the 
homogeneous first order differential equation for dip / dip : 



d f dtp \ 3 dh dip ^ 
dip \ dip J h dip dip 

which has the general solution of the form: 



dip 



1 

c— . 

^ 3 



(A5) 



(A6) 



with c being an arbitrary constant. In order to obtain the 
general solution of the inhomogeneous equation i(A4l) , we 
let c be a function of ip, i.e., dip /dip — c(ip)/h 3 , then plug 
it into Eq. (|A4|I and find 



dc 
dip 



which means 



c{ip) 



= -1, 



-ip + const. 



(A7) 



(A8) 



Combining Eqs. l|A6jl and l|A8j) . we find that Eq. (|A4jl 
has a general solution of the form: 



iP -> C x + C 2 



tp ■ 



a 

r 

(A9) 

where C\ and C*2 are arbitrary constants. We can dis- 
play the divergence of ip> in this limit more explicitly by 
expanding the right-hand side of Eq. \K9i\ in terms of 
(a/2-tp): 



iP 



i(G-f 



dh f a 
dJp V2 



dh /a\ (a \ - 1 

Tp>\2) Kl'V +consi ' 



(A10) 



where we only retained the divergent terms and the con- 
stant term, and accordingly, 



dip 
dip 



dh 

dp 



dh / a\ fa 
~(kp V2/ V 2 

(-- 

\2 



-3 



(All) 



Now it becomes immediately apparent that condition 
(|21|l demands C 2 = a/2, and only the first order di- 
vergence of ip is allowed (which is always present, as is 
clear from Eq. (|A10|l ). 

On reflection, the physical content associated with the 
boundary condition l|21|l may invite further exposition. 
Taking into account Eqs. © and {TQJ) , we note that the 
boundary condition (|21|l states physically hv = at the 
contact line. Consider a region within distance 5 from 
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the boundary. The influx hv should be balanced with 
the evaporation flux, which is proportional to JS. Let 5 
go to 0, and our result follows. Mathematically, one can 
argue that singular behavior of Eq. I|19fl and the order of 
divergence of ip in the limit tp — > a/2 are uniquely de- 
termined by the term 1/h 3 . The solution with a higher- 
order divergence, though compatible with the mathemat- 
ical structure, is not allowed by the physics. 

To make the boundary condition (|21[1 easier to handle 
mathematically, we define a regularized function: 



X = h ip, 



(A12) 



where h 2 is introduced to compensate for the second or- 
der divergence in ip at the contact line allowed by Eq. 
(|19fl . The original problem is converted to the standard 
boundary value problem: 



dh dx 



~ h fx 

dip 2 dp dp 



d 2 h 



2(3v-2)h x=-l, (A13) 



dip 



= 0, 



where x ls defined in the domain —a/2 < p < a/2 and 
is even in ip. The boundary problem (|A13|> has a unique 
solution, which we plot in Fig. 0] for two opening angles 
(q = 7r/4 and Sir/ '4). 

The introduction of the regularized pressure function 
X is not special to the uniform evaporation case, and it 
could be readily defined with a more general and sin- 
gular evaporation profile. Function h in Eq. (| A 12(1 is 
independent of the evaporation, and its exponent 2 is 
introduced to compensate for the possible singularity of 
the reduced pressure function ip in the uniform evapora- 
tion case. With a more general evaporation rate, singular 
behavior of ip may depend on the singularity of the evapo- 
ration rate, and the exponent of h in the definition (|A12|) 
should be adjusted accordingly. It can be shown that in 
the case of the dry-surface evaporation the exponent be- 
comes 5/2 due to the divergence of the evaporation rate 
at the contact line. 

We use both functions ip and x m this paper. Function 
X is employed to obtain the numerical results because 
of its regularity, and function ip is used in asymptotic 
analysis because of its simple asymptotic form and its 
direct connection to the physical properties of the system 
and their singular behavior. 



APPENDIX B: CROSSOVER AT THE 
90-DEGREE OPENING ANGLE 

It has been shown 0,0] that in the limit r <C R(t) w 
Ri the first two leading order terms in the full expansion 
of the surface height h(t, r, ip) read: for the opening angle 



a in the vicinity of tt/2, but not strictly equal to 7r/2: 



h(r, p) 



r ir/a 



1 r 2 / cos2<^ 

7TT ~ -1)+C(a) -, 7 

4 Ri \ cos a / K ' R"/^ 1 <> 



r I nip 
cos — , 



(Bl) 



and at exactly a = tt/2: 

1 r 2 r 
h(r,p) = — In — cos2<^ 

IT Ki Hi 

r 2 (\ 1 \ 
-I \-psm2p hC cos2w , (B2) 

Ri \TT 4 ) 

where C (a) and Co are related by Eq. 1)16[) . and Co is 
independent of a. For simplicity, we replaced R(t) with 
Ri and hence suppressed the time dependence of h. 

As mentioned before, the apparent divergence in Eqs. 
(|14J) . 1)15(1. and hence (|L3 L|) at a — tt/2 is artificial, and, 
as an implicit function of a, the surface shape h(r, p) 
is actually continuous at a = tt/2. This can be readily 
checked by expanding the divergent terms in small pa- 
rameter (tt/2 — a) near a — tt/2 as was done in Refs. 5 
and 17. Moreover, as can be shown by the same method, 
all the derivatives of h with respect to r and p> up to any 
order are also continuous at a — tt/2. 

We can actually estimate the size of the neighborhood 
[tt/2 — Aa,7r/2 + Aa] where the first and the second 
leading order terms on the right side of Eq. 1|B1|) are 
comparable to each other by comparing the dominant 
divergent terms in expressions l|Bljl and l)B2|l . For a < 
tt/2, a = tt/2 — Aa, we have 



1 r 2 cos 2p 
4 Ri cos a 

or, since cos a ~ Aa, 



1 r 2 r 
— —In — cos 2(p, 

TT Ri Ri 



Aa 



In- 



R t 



(B3) 



(B4) 



For a > tt/2, a = tt/2 + Aa, and with C(a) given by Eq. 
ifToJl. one can easily obtain: 



r Tr/a 



1 r 



— ; — -cos2o2~ — In — cos2(p, (B5) 

4Aa^/"-i Y it Ri R, ^' x ' 



or, since 



(r/Rtfl*- 2 « (r/R t y iAa /*, 



Aa 



In- 



Ri 



—4Aa/n 



(B6) 



Since Aa is always small whenever r<ii„ the estimates 
(|B4|) and (|B6|) actually provide compatible results. 

In asymptotic analysis employed in this paper as well 
as in Ref. 6, only the term of the smaller exponent of r 
(it /a ^ 2 according to a ^ tt/2) was retained on the right 
side of Eq. I|B1|) . treating a = tt/2 as the limiting case 
despite a crossover of the exponents at this value. This 
approximation does bring about some subtleties on some 
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occasions, and necessary corrections or modifications of 
our results need to be made. Moreover, although the 
existence of the crossover region [ir/2 — Aa, ir/2 + Aa] is 
uniquely due to the asymptotic form of surface shape h, it 
does manifest itself in other properties of the evaporating 
drop. It is interesting to see how the crossover affects our 
asymptotic analysis of the system. 

In the case of the velocity field, it was found in Eq. 
(|29|l that v v at the contact line would vanish in the limit 
a — > 7r/2. However, intuitively, according to the con- 
servation of the fluid mass, locally the velocity should 
be inversely proportional to the slope of the surface of 
the drop. As mentioned above, dh/dip remains finite 
and depends continuously on the opening angle, even at 



a = 7r/2, and therefore 



V \a=Tr/2 



should not be vanish- 



ing as well. Indeed, if, instead of employing Eqs. i|14|) 
and JTBJ, we use the leading order term h(r,tp)\ a=n / 2 — 
— (1 / n)(r 2 / Ri) \n(r/Ri) cos 2^ on the right-hand side of 
Eq. (IB2|) in Eq. (|2HJ) and let v — 2, we obtain: 



I J a n / r 

UvU=T/2 72 [r~ 



Ri 



(B7) 



which is non-vanishing at the contact line. 

Combining Eqs. (|29J) and (|B7J| . we can again identify a 
crossover region near the opening angle a — 7r/2, which 
was somewhat concealed by the asymptotic analysis we 



employed in this paper as well as in Ref. 6. We believe, as 
partly shown above, that the velocity field (as well as all 
other physical properties of the system) depends on the 
opening angle a continuously. And in the small neigh- 
borhood [tt/2 — Aa, 7r/2 + Aa] actual physical properties 
should be interpolated, so that Eq. H29J) . which holds only 
outside of this neighborhood, is related continuously to 
the result l|B7|l . which applies exactly at a — ir/2. By 
comparing Eqs. (|29fl and (|B7|I . we see that 



tan 



(I 



Aa 











/ 7T 


Ri 



(B8) 



and therefore we can find an estimate for Aa: 



.it 4 
Aa ~ — — arctan — 

2 7T 





7T 










Ri 


~ 4 


Ri 



(B9) 



Result (|B9|I is identical to Eq. i|B4|l . To no surprise, the 
crossover originating from the surface was recovered in 
the result for the velocity, which was determined using 
that surface shape. Similar interpolative estimates can be 
conducted for all other physical quantities. In principle, 
all the results we obtained so far (as well as those in 
Ref. 6) apply only outside of the neighborhood [tt/2 — 
Aa,7r/2 + Aa]. 
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